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Abstract 

a^ . 

' We consider the symmetric simple exclusion process on Z*^, for d > b, and study 

. the regularity of the quasi-stationary measures of the dynamics conditioned on not 

^ I occupying the origin. For each p e]0, 1[, we establish uniqueness of the density of 

+^ ■ quasi-stationary measures in L^((ifp), where Up is the stationary measure of density p. 

. This, in turn, permits us to obtain sharp estimates for PupiT > t), where r is the first 

] time the origin is occupied. 
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1 Introduction 

Let {rjt : t > 0} be the symmetric simple exclusion process on Z'^. In this process, there is at 
most one particle per site (i.e. the state space is VL := {0, 1}^ ), and at rate one the contents 
of neighboring sites are interchanged. The homogeneous Bernoulli product measures, say Up 
with density p G [0, 1], are invariant and reversible for this process. Let r be the first time 
the origin of Z*^ is occupied by a particle. We are interested in two issues: (i) to estimate the 
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probability that the origin remains empty for large time when the initial configurations are 
drawn from Vp for each p g]0, 1[; (ii) to describe the law of 77^, at large time t, conditioned 
on {r > t}, the event that the origin is empty up to time t. 

When the dimension of the lattice is larger than 4, we show that there exists a measure 
/ip, such that for any continuous function / 

limE.^[/(r/,)|r>t]= / fdp^. (1.1) 

This establishes the so-called Yaglom limit [0 . Such limiting measures can be intrinsically 
characterized as fixed points of the semi-groups {Tf,t > 0} defined by 

(r,(/x),/):=E^[/(r/t)|r>t], t > 0. (1.2) 



Thus, fixed points of {T^, t > 0} are dubbed quasi-stationary measures |T^, 0. Here, we 
study the regularity of pp and uniqueness when the dimension > 4. This, in turn, gives us 
sharp asymptotics for the probability of Pup{T > t), namely 

exp(-A(p).t) ^ f(p\\ (1.3^ 

where — A(p) < is the top of the spectrum in L'^ii'p) of L, the generator of the simple 
exclusion process absorbed when hitting {77 : ?7o = 1}. 

We briefiy summarize some relevant results of [Q. In dimensions 1 and 2, the Yaglom 
limit is 5o, the measure concentrated on the configuration with no particle (and A(p) = 0). In 
dimensions 3 and 4, A(p) > for p e]0, 1[, and Ts{yp)ds/t converges to a quasi-stationary 
measure /ip. By analogy with the case of independent random walks [Q, we conjecture that 
the Yaglom limit exists and that /ip is singular with respect to z/p. Thus, it is only for 
dimensions larger than 4 that we expect regularity of /ip with respect to z/p. 



2 Notations and Results 

Henceforth, we consider dimensions larger or equal to 5, and p g]0, 1[. The symmetric simple 
exclusion process (SSEP) on the lattice can be graphically constructed "a la Harris" [[7| 
as follows. First, fix the initial configuration by assigning to each site of l/" a value in {0, 1} 
which indicates if the site is occupied or empty. Then, to each bond -pairs of adjacent 
sites- associate a Poisson processes of intensity 1; Poisson processes of different bonds are 
independent and independent of the initial configuration. At the times events {marks) of 
each Poisson process, the values of the corresponding sites are interchanged. In this way, 
each particle jumps when a mark is present; two particles may jump at the same time in 
opposite directions. By labeling particles, we can trace in time their trajectories: they evolve 
as the so-called stirring particles. This construction is described in Arratia [||. When the 
labels of the stirring particles are disregarded one obtains only the occupation numbers; in 
this case the resulting process, called 77^, has infinitesimal generator 
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where ^'-^ = ■r]k + {5kj - Ski){Vi - Vj) and i ~ j means that |ii - ji| H h \id- jd\ = 1- It is 

well known that the process is Feller, and the product measures of density p in [0, 1], say z/p, 
are reversible for L (see Chapter VIII of Liggett [l^). In other words, L is an unbounded 
self-adjoint operator in L'^[dvp), and local functions form a core for the domain, say D{L). 
We denote by P^^ the law of the SSEP with initial measure Up. Let A = {rj : rjo = 1} and 
denote by r the time of first occurrence of A. As we are interested in the Dirichlet problem 
on A'^, we introduce = {(/?€ L'^{vp) '■ '{'{jf) = for G A}. Let L be the operator 
defined by 

Lf = iAcLf, for feDiL)nnA. 

This corresponds to the simple exclusion dynamics absorbed when hitting the event A. L is 
self-adjoint on Hj, with respect to Up. We call {St,t > 0} the corresponding sub-Markovian 
semi-group of bounded operators on L^(?i_4, Up). In other words, 

Vt>0, SJ{v) = E^[f{r^,)^,^,y]. 

We denote by Ttivp) the probability measure defined by duality on G 7i^, 



where we have used reversibility to obtain the third term. Thus, if ft is the density of Tt{vp) 
with respect to Vp 

It was established in that a non-trivial quasi-stationary measure, say fip, could be obtained 
as limit along a Cesaro subsequence of Tt{vp). Our main result is the following theorem. 

Theorem 2.1 // the dimension is larger or equal to 5, then fip is absolutely continuous with 
respect to Pp. Moreover, for any integer k > 1, dfip/dup G L^{i^p) and 

hm / (^^-^Ydup = 0. (2.2) 



t^oo J y dup dvp 



Remark 2.2 This is stronger than establishing the Yaglom limit, i.e. \im.Tt{vp) = fip. As a 
consequence, / := dpp/dup belongs to D{L) and satisfies (in the L^(z/p)-sense) 

Lf + Xip)f = 0, and SJ = e~^^^^'f, (2.3) 

with (see Theorem 2 of Q) 

A(p) = inf{^^^^^^: feD{L)nnA}. (2.4) 



Theorem |2.1| is based on the apriori bounds through the following lemma, in which we 
reformulate a general result essentially contained in M. 
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Lemma 2.3 Let {St} be the semi-group of a process absorbed when hitting a set A. Assume 
that {St} is reversible with respect to v, and let X < oo be given by ( p.4|) . The following two 
conditions are equivalent 



(i) sup — — < oo, and (ii) sup / f^di> < oo. (2.5) 

t>o "j/(t > t) t>0 



Moreover, if either (i) or (ii) holds, then the Yaglom limit ji exists and ( p.2|) holds. 

Now, the apriori bounds is a corollary of the following proposition, interesting on its own. 

Proposition 2.4 Let the dimension d > 3. Let i G Z^\{0} and rj E Q with rji = 0. We 
denote by rj^ the configuration identical to t] except on i, where its value is 1. There is a 
constant Cd, independent of i and rj such that for any t > 0, 

< P''(r >t)- P"\t >t)< CdP''\T > t)F{Hi < oo), (2.6) 

where Hi denotes the first time a symmetric random walk starting at i hits the origin. 



The relation ( |2.6| ) would be obvious if the particles were independent. Though it is rather 
intuitive for the symmetric exclusion, our proof is rather long. A sketch of it is as follows. 
We first write P^{t > t) in terms of a dual process, say {X(0,i(:)}, which corresponds to a 
stirring process on Z'^\{0} with birth at the nearest neighbors of the origin and with initial 
condition an empty configuration. Then, P^{t > t) — P'^\t > t) corresponds to the weight 
of all paths whose end-points X(0, t) = U U {i} with r]j = for all j eU C Z'^\{0, i}. The 
problem is then to uncouple U from {i}. We then re-express P(X(0,t) = U U {i}) in terms 
of a dual with finitely many particles, say {A^}. Note that {A^} is not the 'natural dual' 
of {X(0,t)} and the correspondence is obtained through a Feynman-Kac formula. Then we 
show a correlation inequality for the expression in terms of {A^} by generalizing Andjel's 
inequahty 

The results about apriori bounds is the following. 

Corollary 2.5 Let the dimension d > 5. (i) There is a product measure of density a{i) 
for i such that for any t > 

-< Tt{vp) -< z/p, and ^^(1 - — f < oo, 

where -< denotes stochastic domination. 

(ii) For any integer k > 1, there is a positive constant C , such that 

sup / fhv)diyM < C. (2.7) 
A consequence of Lemma |2]^ is a sharp asymptotic estimate for the tail of r (compare with 



Lemma 1). 
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Corollary 2.6 // the dimension d > 5, then 

„-\(p)t r 

lim — = / fdvp. (2.8) 

Finally, we have a uniqueness result and some properties of fip. 

Theorem 2.7 (i) The set {/i <^ z/p : fi is quasi stationary and dfi/dup G L^(z/p} contains 
only fip. (a) For z G Z"', define 9i : Q ^ Q with Oirjk = rjk+i, then for any (p local (i.e. 
depending on finitely many sites) 

lim / (p{eir])dfip{r]) = / cpdup. 

ll'lHooJ J 

(Hi) If u is a probability with a continuous density with respect to Vp, then 

1 /•* 

hm - / Ts{y)ds = fip. 

t^OO t Jq 

The convergence holds in weak-L?'{vp). 



Proposition p^.4| is proven in Section ^. Section ^ contains the proofs of Corollary |2.5| , and of 



Theorem |2.1| . In section ^, we establish the uniqueness part of Theorem In section ^ 
we show that in Hp the density at infinity is p, and we conclude with the result about the 
basin of attraction of /Xp. 



3 Proof of Proposition |2t^. 



3.1 Duality and Feynman-Kac. 

We first express P''(r > t) using the dual process ([[l^], 01) based on the fact that the Poisson 
clocks associated to bonds are invariant by time reflections. The dual process tracing back- 
in-time the positions of the stirring particles touching the origin can be described using 
the graphical construction at the beginning of Section 2. Again at each bond, there is an 
independent mark process corresponding to the realization of a Poisson process of intensity 1. 
At each mark between and one of its nearest neighbor, say i, a particle is born at i unless 
i is already occupied (in which case nothing happens); the particles born in this way evolve 
afterwards as stirring particles on Z'^\{0}. The only difference with the previous construction 
is that now it is imposed to the origin to be occupied at all times — so that when it becomes 
empty, it is immediately occupied with a newly created particle. Assume that at time 0, the 
lattice is empty and let X(0,t) be the set of sites occupied by the stirring particles at time 
t; all these particles have been created at the origin. Let P denote averages over the Poisson 
realizations. If V* is the collection of finite subsets of Z'^\{0}, then the duality formula reads 
for any t > 

P^{r>t) = (1-r^o) E^(^(^'^) = ^)n(l-^^)- (3.1) 

AeV* jeA 



5 



Thus, if rj is such that rji = 

P\r>t)-P^'{T>t) = (l-vo) ^{X{(ll,t)=A)l[{l-Vj). 

Assume, for a moment, that for U eV* and i ^ U, we have a constant Cd independent of i 
and t such that 

P (X(0, t) = UU {i}) < CdP (X(0, t) = U) F{Hi < oo), (3.2) 

where we denote by Hi the first time a symmetric random walk starting at i hits the origin. 
Then, for 7] such that rji = 

P\T>t)-P^\r>t) < CdnH^<oo){l-r^o) Yl P(^(0,^) = f^) H ' V,) 

uev*,i^u jeuu{i} 

<Cd{l~r]i)¥{Hi<oo)P'^\T>t). 

Thus, it remains to prove (^.21). 

Let be the generator of {X(0, t), t > 0}, and let be the associated semi-group. We 
first express the dual of {X(0,t), t > 0} in terms of a process with finitely many particles. 
Actually, we are only interested in 5'i+(lA)(0) := P(X(0,t) = A) for A G V*. Let A and A 
be in V*. We have, using A for the symmetric difference, 

^^{Ia){A) = Yl lUAA{x, y}) - U{A)] + U {y}) - 1a{A)]. 

x^y; x,y^{0} ?;~0 
\AA{x,y}\ = \A\ 

The first sum corresponds to the stirring process over Z'^\{0}, while the second sum corre- 
sponds to birth at the origin. We reexpress now the last sum. For simplicity, we omit to 
write ?/ ~ 0. Thus, 

E[1a(^ U {y}) - U{A)] = YM^ U {y}) - UA)] - Y U^)- 

y^A y^A y^A 

t/SA y^A 

We claim that this expression is equal to 

C := 5^[U(A\{i/}) - U(A)] - Y + Y 

j/SA j/^A yeA 

Indeed, we expand C 

C = E[l^(M{l/})-lA(A)] + $^[U(A\{y})-U(A)]-5^U(A) + J]U(A) 

yeA j/eA y^A yeA 

y^A yeA 

= E[i^(^\M)-1a(A)]-5^U(a). 

y<^A y^A 
yeA 
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Thus, calling Ao := {y G Z'^ : y ~ 0} and using the self-duality of the stirring part of 
we obtain 

^+(1a)(^) = E [1a(AA{x,i/}) - U(A)] + J2l^A{A\{y}) - U{A)] + V{A)1a{A). 

\AA{x,y}\ = \A\ y&A 

(3.3) 

where we set V{A) := 2|A/o fl A| — |A/o|. Now, let C denote the generator of the stirring 
process on 'Z'^\{0} with death when particles jump on the origin. Then, ( p.3|) can be written 
like 

C+{Ia){A) = C-{1a){A) + V{A)1a{A). 
Thus, if u{A,t) := 5+(1a)(0), we have 

^^^^ = 5+(/:+1a)(0) = C-u{A,t) + V{A)u{A,t). 
Now, we call {A{t),t > 0} the process generated by and we use Feynman-Kac to obtain 

P(X(0,t) = A) = E^[e-fo^^^(')^''n^{A{t))] = e-l-^«'*E^[exp(2 / \Afo n A{s)\ds)l^{A{t))]. 







Let now U E V* and i ^U. We show in section ^ that if A G P* and 



^(A,t):=E^[exp(2 / |AronA(s)|ds)l0(A(t))], then giUU{t},t) < giU,t)gi{t},t). (3.4) 



Then, in section ^]3| we prove that g{{i},t) < CdP{Hi < oo) for d > 3. Inequality 
follows then readily. 

3.2 A generalized correlation inequality. 

To make the notations closer to those of Andjel [|I|, we set p{x,y) = 1 when x ~ ?/, and 
p{x, y) = otherwise. Also, we realize our stirring process as an exclusion process: the 
particles attempt to jump to one of their nearest neighboring sites at the time marks of 
independent Poisson processes of intensity 2d; if the site chosen (each neighboring site is 
chosen with the same probability) for the attempt is occupied, the particle stays still. As we 
are blind to the labeling of particles, the trajectories are, in law, indistinguishable from our 
initial stirring process. 

We proceed by induction on n to prove that for any sets A, B E V* , with An B = ^, for 
any t > 0, a G M, and any n-tuples < < < ■ ■ ■ < s" < t, the following inequality 
holds 

(3.5) 

This will be our induction hypothesis at order n. Once (|375|) is proved, inequality (|3T 
follows easily as in Proposition 4.1 of M. 



7 



step n=0. We need to prove that for A,B eV* with A n B = 

P^uB(^(^) ^ 0) < pA(^(^) ^ 0)P^(A(t) = 0). (3.6) 

Following an idea of Arratia we represent the process K{t) as limit of a stirring process 
with no absorption on an enlarged lattice: we link the origin with 0, the origin of a 
three dimensional lattice 1? isomorphic to 1? (here we fix for concreteness; any graph 
supporting the stirring construction, for which the corresponding random walk is transient 
would fit). On each bond of 1? and on the bond (0, 0), the rates of stirring are set equal to n 
large. On the enlarged lattice 7/" U Z^, the particles perform a conservative stirring process, 
though with different rates whether they jump across the bonds of l/' or across the bonds 
of 7? and (0, 0). When a particle hits the origin 0, it has a probability going to 1 as k — *• oo 
to wander in 7? up to time t without using bonds of Z^. We call U{t) the stirring process 
on Z'^ U Z^, and the law of the Poisson marks on the enlarged lattice. It is not difficult 
to show that for any A G P*, 

(A(t) = 0) = lim (t/(t) C I?). (3.7) 

Now, for the stirring process on the enlarged lattice Z'^ U Z^, we use a correlation inequality 
due to Andjel [0]: 

P;f^^(f/(t) C I?) < ^i{U{t) C I?)¥^{U{t) C I?). (3.8) 
Thus, after taking the limit k oo, we obtain ( p.6|) . 

Step n. Our proof follows essentially Andjel's proof. Our induction hypothesis is that ( |3.5| ) 
is valid for n — 1 instants of time. Let < < < • ■ ■ < s" < t be ra instants of time, and 
for each A eV* let 

n 

g^A, t;s\..., O = E^[exp(a {Afo n A{s%mA{t))]. 

i=l 

We set A = 2(i(|y4| + |-B|) and we let ri be the first time a particle of AU B attempts a jump 
(i.e. Ti is an exponential time of parameter A). Note that by the Markov property 

n 

E^^^[exp(a J2 1-^0 n A{s%U{A{t))^,^^,^}] 

i=l 

n 

= Pin > s^)e"(l^'^«l+'^"^°l)E^^^[exp(«^ lATo n A{s' - s')\)mA{t - s^))] 

i=2 

= P{n > si)e"(l^^°l+l^^^°l)(7„_i(A UB,t-s';s^-s\..., - s'). 

Following III, using the shorthand notation s for s^, . . . (and its abuse s — u = (s^ — 
. . . , — u)), and writing AA{x,y} to mean AA{x,y}\{0}, for we have to account for 
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deaths of particules when they jump on 0, we have 

gn{A U 5, t- s) = P(ri > si)e°l^"^ole"l^"^«l(?„_i(A U 5, t - - s\ . . . , s" - s^) 

+ [ y)9niAA{x, y}UB,t-u;s- u)] 

x&A 
y^AuB 

+ [ p{x,y)gn{AU BA{x,y},t - u;s - u)] 

y^AuB 

+ [Yp{x,y) + Yp{x,y)]gn{AUB,t-u;s-u)^ 

xGA xeB 
y&B ydA 



Reasoning as if the particles in A were independent from the particles in 5, we obtain 

gn{A,t] s)gn{B,t]s) = 

P{n > sOe°l^"^°le'^l^^«l(?„_i(A, t - s'; - s\ . . . )gn-i{B, t - s'; - s\ . . .) 

+ / rfwAe"^"- 1 [ ^ p{x, y)+ y)]9n{A, t-u;s- u)gn{B, t-u;s-u) 

x,yeA x,y£B 

+ [ X] y)9n{.AK{x, y},t-u;s- u)gn{B, t-u;s-u)] 

xGA 
y^AuB 

+ [ X] t-u;s- u)gn{BA{x, y},t-u;s- u)] 

y^AuB 

+ E ?/)5'n(^A{x, y},t-u;s- u)gn{B, t-u;s-u)] 

xGA 
yeB 

+ l^Pix, y)gniA, t-u;s- u)^„(5A{a;, y},t-u;s- u)]|. 

xGB 

yeA 



Define 

Gn{t) = sup sup 5(„(CUL',t;s) -5(„(C,t;s)5(„(L',t;s), 

0<si<---<s"<t cnD=0 

also set = sup{G„(s) : < s < t}. Now, the key observation of Andjel, in p. 720, 

is that for x & A and y E B (so that both x,y ^ 

gn{A U B, t; s) < F„(t) + ^ t; s)^„(5A{x, y}, t; s) + ^„(AA{x, y}, t; s)^„(5, t; s)) . 
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Thus, using the induction hypothesis = 0) and the symmetry of .), we obtain 

gn{A U B, t- s) - gn{A, t- s)gn{B, t; s) 



< e 



-\u 



^p{x, y) y) +^ p{x, y) +^ p{x, y) + 2 ^p{x, y) 



x,yeA 



x,yeB 



xeA 

y^AuB 



xeB 

y^AuB 



xeA 

yeB 



Fnit — u)du 



= f e-^''XFn{t-u)du<Fn{t) [ Xe'^'^du. 
Jo Jo 



Thus, by taking the supremum over A,B e V* with A (1 B — we obtain Gn{t) < 
Fn{t) Xexp{—Xu)du. Finally, this imphes that F„(i) = 0, and the proof is completed. 



3.3 Upper bound gi{i},t) < Cd^Hi < oo). 

We first use the classical representation of the trajectories {A{{i},t),t > 0} as sequences 
of jump times {Ti,i e N}, which are independent exponential variables of parameter 2d, 
associated with the paths of a simple symmetric random walk killed at the origin, say {Aj, i e 
N}. The processes {Aj, i e N} and {r,-, j G N} are independent. We use the notation and 
to denote average over paths of {Aj, i G N} starting on y. Let Tq = inf {n > : A„ G Ao} 
with A/q = {y : y ^ 0}. When Tq < oo, let Ti be the first return time to A/q, whereas when 
To = oo, set Ti = oo. Then, the sequence of successive entrance times in A/q, {T2, T3, . . . } is 
defined inductively. The number of return to A/q is called R, i.e. if Tj = 00 but Tj-i < 00 
then R = i. Our walk on being transient, we have R < 00, a.s. We note also, that by 
symmetry, for any y, y' G Ao 

P^(ro < 00) = P^'(ro < 00) and VA;gN, py{R ^ k) ^ P^' {R ^ k). 

For convenience, we call P^o(To < 00) := py{% < 00) and P^°{R = k) = py{R = k) for 
y eMo. Now, 

R 00 k 

g{{i}, t) < [l{ro<oo} exp(2 ^ ttJ] =J2^' [l{ro<oo}l{K=fc} exp(2 ^ rrj 

i=0 fe=0 i=0 

00 

- P\To< 00) J2P^"iR = k){E[e^^'])\ (3.9) 

k=0 

where we used the Markov property and induction. Now, by the same arguments 

P^°iR ^k)< (P^°(To < 00))'' . 
On the other hand, the evaluation of £^[exp(2Ti)] is easy 

E[e^-^] = r e^'2de-''^'dt = (3.10) 
Jo d — 1 
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Thus, with ( |3.9| ) and ( p.lO| ), our upper bound follows easily as soon as 

P^o(To < oo) < (3.11) 

We want to formulate ( p.ll| ) in terms of hitting probabilities for the standard random walk, 
say {Sn,n > 0}. We will denote the averages over the standard walk with a tilde. Let 
K = inf{n > : Sn = 0}, and note first that 

P^"(To < oo) = P^'X^o < oo, K > To). 

By conditionning on the first move, we obtain 



P-^«(To < oo) = ^ + P^»(To < oo, K > To 



2^ , ^ . oj 

Thus, ( pUP is equivalent to P-^o(To < oo) < {2d - l)/2d. We recall that R is the number 
of return to A/q for a walk starting on Ao- We can set 5*0 = but count what happens only 
after two steps 

oo 

^ = X] ^{-^ne-^o} = X] l{r„<oo}- (3.12) 

n=2 n>l 

So that 

ti 1-F«.(r„<a)) 
Thus, ( p.ll|) reads E[R] < 2d - I. Now, we note that for n > 

P{S„.+i = 0) = P(S„,+i = 0,S^ = y,ye K) 

= J2PiSn+i = 0\Sn = y)P{Sn = y) = j^P{SneAfo). (3.13) 

Thus, 

oo oo 

J]P(5„ e Afo) = 2dJ2P{Sn = 0). 

71=2 n=3 

In dimension 3, it has been established (see page 170, exercise 2.7) that 

oo 

1 + J2 Pi^n = 0) = (V6/327r3)r(l/24)r(5/24)r(6/24)r(ll/24) = 1.516... 

n=l 

Thus, in dimension 3, E[R] < 6(0.52) < 5 and our condition ( p.ll| ) holds. Thus, in dimension 
3, E[R] < 6(0.52) < 5 and our condition ( ^.11[ ) holds. We conclude that ( |3.11| ) holds for any 
dimension larger or equal to 3 because the right hand side of E[R] <2d — l increases and the 
average number of visits to decreases with dimensions. As we have not found a reference 
of this latter fact, we present a short proof due to Andjel. The number of visits to is a 
geometric random variable of parameter Po(To < oo), so it suffices to show monotonicity 
for this quantity. Let i be any neighbor of 0; by symmetry a{d) := Po(Tq < oo) = Pi(To < 
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oo), for dimension d > 3. Project the (i-dimensional walk on a hyperplane orthogonal to 
i', a neighbor of the origin different from i. The projected walk on Z"^'^ has transition 
probabilities l/2d to go to each of its 2{d — 1) neighbors and 1/d not to move. It is clear 
that a{d) is not larger than the probability of visiting the origin starting at (the projection 
of) i for the projected walk. This latter probability is actually equal to a{d — 1). Indeed, 
the projected process goes along the same trajectories as the {d — l)-dimensional standard 
walk and the waiting times at each point are geometric random variables with parameter 
[d — l)/d: thus, if the trajectory of the standard walk is such that {T < oo}, then the same 
holds for the projected walk and vice-versa. 

4 Apriori Bounds. 
4.1 Proof of Corollary [275| 

(i) The proof proceeds along the same lines as the proof of Theorem 3c) of once the 
measure z/q,( ) is defined. We set = 0, and for ^ 7^ 0, let be defined by 

ftj 1 — p 1 

P ^ l + CdFiHi < 00)' 

where the constant Cd is that of ( |2.6| ). Note that < < p. Now, in the proof of Theorem 
3c) of @], we showed that z/a(.) <^ z/p and dh'a{.)/di'p G L^^Up) for p > 1 as soon as 

— — )^ < 00, or equivalently in d > 3 ^{Hi < 00)^ < 00, 

which holds as soon as c? > 5. 

We rewrite ( p.6|) on {rj : rji = 0} with i ^ 0, denoting by CTj the action of spin flip at site 
i (i.e. {crir])k = r]k if k ^ i and {air])i = 1 - r]i) 

adt>j^^ft or equivalently a,(/,)^ > a,(^)/,. (4.1) 
1 — ctj p dUp dUp 

Now, on A, we form ip = dTt{h'p)/di'a{.) and we note that ip is increasing. Indeed, if i 7^ 
and rji = 0, then i ^A\j is nothing but ai(p > ip. Now, as a product measure i'a(.) satisfies 
FKG. Thus, for ip increasing 

j ipdTtiup) = j i)(^^^)dva{.) > J tpdua(.) j = J 

Thus, z/q(.) -< Tt{up). The fact that Ttivp) -< Up comes from the fact that ft is decreasing 
and Vp satisfies FKG. Now, (ii) of Corollary p. 5| follows as in the proof of Theorem 3c) of [0]. 
Using that ft and dh'a(.)/di'p are decreasing, for 2 > 1 and j > 0, 

P' - P'-' (t^) <- (^) ^ P'- i^T"- 
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Thus, we obtain by induction, for each n > 1 



fl^du, < / ( ^ ) du,. (4.2) 



'dUa{^ 

Since the right hand side of ( |4.2|) is bounded for d > 5, the corollary follows. 
4.2 Proof of Lemma Li.'6 



In 0] section 4, we have that t ^-^ R{t) := e / Pyij > t) is increasing and -R(O) > 0. 
Suppose (i) and let limt^oo -R(''^) = R < oo. Now, (ii) follows from 

_ f St{lA^)St{lA^) , _ f S,,{U^) _ P^jr > 2t) _ R(2t) R 

J PAr>ty J PAr>tr'^~ PAr>ty~ Rity- RW ^ ^ ^ 

Conversely, Let /i be a limit point of {1/t Ts{h')ds} along the subsequence in weak- 
L'^{dv). By Theorem 1 of [Q, /i is a quasi-stationary measure with P^{t > t) = exp(— At). 
Thus, {€} follows from 



lim Rit) = lim — / — — -ds = lim — / / fs- — dupds = / - — dup < oo. 

t^oc n^+ootnJo Pvp[r > s) n^+oo J Q J dUp J \dlJpJ 

(4.4) 

Now, if either [i) or {ii) holds, by remark 3 of 0], we have that the Yaglom limit exists. 
Now, to show that j {ft — f ^dv converges to 0, we only need to show that j ffdv converges 
to / pdv. It is easy to see from (g]3[) and (^), that lim J f^dv = limi?(t) = J {dn/duYdu, 



and the proof of Lemma 2.3 is concluded 



5 Proof of Theorem |277| (i): Uniqueness 
5.1 Positivity of /. 

We first show that Up-a.s., / > on {?7o = 0}. We introduce a symmetric simple exclusion 
process on Z'^\{0}: there is no site 0, and its adjacent bonds are suppressed. Let v* be 
the product Bernoulli measure on Z'^\{0} of density p. Let {5;, t > 0} be the Markov 
semi-group of this process. It is known that the process {S^, t > 0}, with initial measure 
I'*, is reversible and ergodic: indeed, by Theorem 1.44 on page 377 of |]TD|], i^p is an extremal 



invariant measure and by Theorem B.52, on page 23 of [[T^], all extremal invariant measures 
are ergodic. In other words, if for any t > 0, S^g = g, i^*-a.s., then g is constant t'*-a.s. 

Let B := {t] : TjQ = 0, /(r^) = 0}, and note that Up-a.s. 

Stf.lB = exp{-Xt)f.ls = 0. 
Now, because / > 0, we have that for any e > 0, el{/>e} < /• Thus, for t] & B, Up-a.s. 
We > 0, P^'iifivt) >e},T>t)=0^ P^iifiVt) > 0}, r > t) = 0. 
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Now, if the bonds hnking to its neighbors are not marked, up to time t, then St acts hke 
S^. So, if To is the first time a Poisson mark appears in one of these bonds, for rj e 8, Up-a,.s. 

P^{f{Vt)>0,To>t) = Q, 

because Tq < r. Now, Tq is independent from the Poisson processes on bonds of Z'^\{0}, so 

P'^ifiVt) > 0, TO > t) = 5;i{;(.)>o}P(ro > t) = 0. 

Now, for any t > 0, P(to > t) > 0. Thus, for rj e B, u*-a.s. 

StMfi.)>o}{v) = 0. 

In other words, for any i > 0, we have v*-a..s. 

s:ib > Is- 

Now, u* is reversible for S^, so both expressions have the same mean, and we conclude that 
i/*-a.s., for any t > 0, S^Ib = Ib- By the ergodicity of u*, we conclude that 1b is i^* constant, 
so that necessarily Up{B) — 0. 

5.2 One eigenvalue with a positive eigenvector. 

Suppose that /, /' G 7i_4 are the densities of two quasi-stationary measures. There are two 
real numbers A(p) and A' such that /, /' satisfy in an L'^{i'p) sense 

^t/ = e-^^'')*/, and ^,/' = e-^V- (5.1) 

First, we show that A(p) = A'. We assume that = fdvp corresponds to the Yaglom limit. 
Thus, 

^hm j fftdup = J f'fdup, and ^hm J f f^dup = J fdup. (5.2) 
However, as / > on {770 = 0}, i/p-a.s., 

e-Xr^ . /ii^ . / ff,,.^ _ / ff,.^ > 0. (5.3) 
P,(t > t) P,(t >t) j J J' P J •' ' ^ ^ 



Similarly, 



Thus, A(p) = A'. 



/ 



-X(p)t 



jfdi 
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5.3 Dual expansion. 

We expand / on the countable basis of L'^{vp), say {Ha-, A G P}, where V is the collection 
of finite subsets of 7/- and 

H, = l, and HA{ri) = \{-^^^=^. (5.5) 
Thus, there are real numbers {Ca-, A G V} with 



f=J2 CaHa, and f fdvp =Y.C, 



(5.6) 



The constraint that / G Ti^^ i.e. rj^firf) = 0, is equivalent to 



VA ^ 0, J HAVofdu, = 0. (5.7) 

Thus, 

Vp(1 - P) y HAumfdu, + pj HAfdu, = 0. (5.8) 

So, for A ^ 0, 

v/p(l - p)Cau{0} + pCa = 0. (5.9) 

We define for all A e V* 

= (-y^) Ca. (5.10) 



Now, condition ( |5.9|) reads '?/'(AU {0}) = ip{A) for A ^ 0. Now, we express {—Lf, /) in terms 
of the {Ca}- 

{-Lf,f) = J2J2^Cs-CAf, (5.11) 

Aer B^A 

where B ^ A ii there is i G A\B and j G B\A with AAB = {i,j}- We replace the Ca's by 
the '?/^(yl)'s and distinguish to ehminate ( |5.9| ) 



-^/'/) = E( E 7l^l(^(i?)-^(^))'+ E 7l''l(^(i?\{0})-^(A))M 

A^O \B^O,B'^A B30,B^A J 

+ E( E 7l^l(^(5\{0})-V^(A\{0}))^+ 5^ 7l^l(^(5)-^(^ 



5^ ^1^1(1 + ^)(V;(i?)-V;(A))2+ 5^ ^\^\{^{B)-ij{A)y 
A^o BeAfA BeMx 

+ E l^^^~'\m)-i;iA)r], (5.12) 



where 7 = (1 — p)/p, S G A/'^i means B A and B ^ 0, B & Ma means 5 ~ A U {0} and 
5^0, and 5 G means S U {0} ~ A and 5^0. With this rewriting, {-Lf, f) can be 
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thought of as the Dirichlet form, £{'ip,%p), of a dynamics with finitely many particles, with 
creation and annihilation at site 0, with respect to a measure m(A) = 71"^!. The advantages 
of this rewriting are threefold: (i) the constraint ( p.9|) has vanished, (ii) the new dynamics 
is clearly irreducible and (iii) in studying the minimizers of S{ip,ip), we can assume the 
{il^i^A)} to be nonnegatives. Indeed, note that £{ip,ip) > S{\ip\, {ipD and equality holds only 
if ip{A) > for each A E V*. Also, we rewrite the LF'{vp) norm of / in terms of the 
{^{A),AeV*} 

\f? = Y.^l + clum = (1 + 7) = (1 + (5.13) 

A^O A^O 

Now, 1/ + A/ = implies that for A e V* 

(^(B) - ^{A)) + 7 5^ im) - m)) + E (^(^) - ^(^)) = -^^i^)- (5-14) 

BeAfA BeM+ BeMx 

Thus, ijj{A) > for al l A e V*. Now , let (f) = ip'^ and note that for any A and B E V*, 
the functional ^-^ {\J (j){B) — (j){A)Y is convex. Assume that and ip' are two positive 
normalized minimizers, and let (f) and (j)' be their respective squares. Then, for any A G [0, 1] 

V'A := VA0 + (1 - A)0' has I IV^aIU = 1 and £(V^a, V^a) < A^(V^, V^) + (1 - A)^(V^', V^')- 

Thus, the convex inequality is an equality, so that for any A eV* and any B E A/'^UA/'^ UA/"^ 

(MB) - = mB) - + (1 - \){fiB) - 

which implies, after expanding, that (f){A)(j)'{B) = 0(i?)0'(y4). Now, as (f){A) > and the 
dynamics is irreducible, we conclude that 4> = 4>' and so are the positive square-roots ip = 

6 Proof of Theorem |277| (ii) and (iii) 

6.1 Proof of Theorem |2.7| (ii): Density at infinity 

The facts that for any t > 0, i'a{.) -< Tt{vp) -< Vp with a{i) p when ||i|| — oo implies that 
for any A E V 

for rijGA ^j+j increasing function. Now any local function cp can be written as a linear 
combination of local increasing functions, and the property follows by linearity. 

6.2 Proof of Theorem |2.7| (iii): Basin of Attraction 

We show that for any measure z/, any subsequence of the Cesaro limit of {Tj(z/)} contains a 
further subsequence converging to a quasi-stationary measure, say /i. When the density of p 
with respect to Vp, say 0, is continuous, we show that p = fip. As in the proof of existence of 
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a quasi- stationary measure (see [^] Lemmal), we establish first the existence, for any s > 0, 
of the following limit 

I^ ^^Pt'^lV^ =^M-Kp)s). (6.1) 
t^oo rv(r > t) 

Indeed, recall that 

Pu{r >i) = j St{lA^)lA4du^ = j St{lA^<P)dvp, (6.2) 
so that by the existence of the Yaglom limit 

[r > t) J 

Thus 
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hm ^5i(L^l±£) = lin, "^--(7^+;^ = exp(-A(p).). (6.4) 

By the weak* compactness, for any sequence {tn}, there is a further subsequence (still named 
{tn} for convenience), and a probability measure /i, such that for any continuous function ip 

^ j\Tt{u),^)dt-^ j ^dfi. (6.5) 

The same argument as in the proof of Theorem 1 of ^ implies that /i is quasi-stationary. 
Now, we establish apriori estimates 



dv, ) J \ P,{T>t) J ' 



This quantity is bounded by Corollary p.5| and ( [6.31 ). By standard arguments, this implies 
that /i G L'^{i'p) which by the uniqueness result establish that /x = /Xp, so that the Cesaro 
limit exists and is /ip. 
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